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-»«. ■ Abstract 

^. 

Q^i, We present an overview of the research activities of the theoretical particle physics 

Q_ij group at the University of Mainz on the calculation of massive one- and two-loop 

Feynman diagrams. The main objective of this research was to develop an automatic 
one- and two-loop calculation program package. A first version of such a program 
was recently realized by Briicher, Franzkowski and Kreimer. We describe in some 
^ ■ detail the present features of this automatic loop calculation program as well as the 

integration techniques that go into the program. The program runs under the name 
of XLOOPS. The present version is labelled by XLOOPS 1.0. The program allows 
one to calculate massive one- and two-loop integrals in the Standard Model including 
their tensor structure - all at the click of a mouse. Renormalization features are in- 
cluded in that the UV divergences in UV divergent integrals are explicitly computed 
in dimensional regularization. One-loop integrals are calculated analytically in d ^ 4 
dimensions whereas two-loop integrals are reduced to two-fold integral representa- 
tions which the program evaluates numerically. We attempt to provide a synopsis 
of the novel loop integration techniques that have been developed for the XLOOPS 
program. They allow for a fast and efficient evaluation of massive one- and two-loop 
integrals. We discuss Higgs decay at the two- loop level as a first application of the 
novel integration techniques that are incorporated into XLOOPS . 



1 Introduction 



In the last few years a large part of the theoretical particle physics group at the University 
of Mainz has been involved in the computation of massive one- and two-loop Feynman 
integrals. This activity was part of a continuing group seminar series on the evaluation 
of Feynman diagrams. Apart from the work of the senior members of the group many 
beautiful new results have been achieved by our diploma and PhD students in the context of 
their thesis work. Present members of the group are L. Briicher, M. Fischer, J. Franzkowski, 
A. Frink, S. Groote, V. Kleinschmidt, J.G. Korner, R. Kreckel, D. Kreimer, A.J. Leyva, M. 
Mauser, K. Schilcher and J.B. Tausk. Past members of the group include A. Czarnecki, U. 
Kilian, K. Melnikov, D. Pirjol, P. Post and O. Yakovlev. We have greatly benefitted from 
the many visitors and guest scientists who have visited our group over the years and who 
have generously shared with us their insights into the calculational techniques of Feynman 
diagram computations. They are too numerous to be all quoted here in name. All of 
this effort would not have been possible without the support of various funding agencies. 
These include the DFG (Deutsche Forschungsgemeinschaft) with individual grants as well 
as stipends provided by the Graduiertenkolleg "Teilchenphysik bei hohen und mittleren 
Energien", the BMBF (Bundesministerium fur Bildung und Forschung) and a HUCAM 
project granted by the European Union. 

This write-up is certainly in the character of a workshop write-up. We take this op- 
portunity to present an overview of the activities of our group and the results obtained 
therein. In particular this means that we do not attempt to give a balanced account of 
the present state of art of loop calculations the world over. We ask for forebearness. It is 
clear that the results that we are presenting are too numerous to have been obtained by 
any single person alone. The results represent the efforts of various subsets of the members 
of the group. When discussing particular results we shall not always specify the subset 
of members of the group that were involved in a particular calculation in order to avoid 
excessive cross-referencing. The respective members can be identified from the references. 

By means of introduction we want to remark that there are two main approaches to 
doing loop calculations. The one approach we would like to call the absorptive-dispersive 
approach. One first calculates the absorptive parts of a Feynman diagram and then uses 
dispersion relations with possible subtractions to obtain the dispersive part of the Feynman 
diagram. The experience is that the absorptive-dispersive approach usually goes further 
when calculating two-loop diagrams if the aim is to perform analytically as many of the 
necessary integrations as possible. This is particularly true for two-loop diagrams with 
internal loop masses. 

There are many examples of two-loop calculations that have been done using the 
absorptive-dispersive integration technique. For example, in Mainz we have calculated 
the two-loop gauge boson contributions to the decay H — > 7 + 7 using the absorptive- 
dispersive technique |T[. The calculations were done using the Equivalence Theorem. In 
the high energy limit the transition H — > 7 + 7 is dominated by the longitudinal degrees 
of freedom of the internal loop gauge bosons which can be represented by the massless 
and spinless Goldstone excitations of an effective scalar theory. The Higgs boson is the 
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only massive particle in the theory. It appears both as external and internal loop particle. 
There are a number of two-loop diagrams that contribute to the process. They are too 
numerous to be all shown here. All absorptive contributions were determined analytically, 
and all but one of the remaining dispersive integrals were done analytically. In Fig.[l| we 



Figure 1: Three-particle cut of a crossed ladder diagram contributing to H — > 7 + 7. Solid 
lines: mass Higgs particles; dashed lines: massless gauge bosons of the effective theory. 

show the contribution of the three-particle cut of the the massive two-loop crossed ladder 
three-point function which was the only contribution which did not yield to analytical 
integration in the dispersion integral. The remaining one-dimensional dispersion integral 
of this one particular contribution was then done numerically. One of the outcomes of 
the analysis was that the results of the effective theory should be reliable in the range 
0.6 TeV < niu < 1.5 TeV. Above 1.5 TeV the two-loop result blows up. 

In another example the absorptive-dispersive approach was used to calculate the two- 
loop contributions to the K° charge radius in chiral perturbation theory 0. Quite naturally 
the loop masses in chiral perturbation theory cannot be neglected since chiral perturbation 
theory is an effective low energy theory. The final result of the calculation was obtained 
in terms of a one-dimensional integral representation which was evaluated numerically. 

As a curious side remark we want to emphasize that it is even sometimes advantageous 
to take the opposite route and first calculate the full loop contribution directly and then 
extract its imaginary part. For example, rate calculations can be done quite efficiently 
in this way. If one shies away from phase-space integrations but knows how to do loop 
calculations e.g. by Feynman parameter methods then one can obtain rates very efficiently 
from the imaginary part of the full loop contribution. This technique has been used to 
great advantage in the calculation of the inclusive semileptonic differential rate of polarized 
Afe's to polarized r's within the Operator Expansion method in HQET ||. The relevant 
one-loop diagram is shown in Fig.^| together with its absorptive cut which determines the 
rate of the decay. 

A large part of the remaining discussion in this presentation will concern itself with 
the direct computation of massive two-loop integrals not using the absorptive-dispersive 
approach. As a general feature of the direct approach we shall find that all but two of 
the massive two-loop integrations can be done analytically. In the absorptive-dispersive 
approach there are many examples where one can go further and where all but one of 
the massive two-loop calculations can be done analytically. A noteworthy example is the 
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Figure 2: Rate of A[ — > X c + r ^ + v T extracted from the absorptive part of the one-loop 
amplitude r +(T) + aJ, T) -> + A[ T) . 

efficient one-dimensional integral representation written down by Bauberger and Bohm for 
the general massive two-loop two-point function M. This has to be compared with the 
two-fold integral representation derived in the direct approach to be discussed later on. 
However, the absorptive-dispersive approach involves a detailed consideration of the two- 
and three-particle cut structure of the contributing diagrams including possible subtraction 
terms and anomalous threshold contributions which is not always simple theoretically. In 
particular the requisite analysis does not possess the flexibility which is needed for a general 
purpose program package that we were aiming for. 

2 Direct Approach to Loop Calculations 

In Mainz we have developed a completely new approach to the calculation of massive one- 
and two-loop diagrams. These radically new integration techniques have been devised in 
a series of ingenious papers by Kreimer || ||, ^ ||. In order to appreciate the quality of 
innovation of these new techniques let us briefly recapitulate the main features of the more 
traditional approach to loop integrations. The essential integration features employed by 
most practitioners in the field may be summarized by the following statements 

• integrate covariantly 

• introduce a set of Feynman parameters 

• perform a Wick rotation and integrate in Euclidean space 

• do tensor integrals a la Passarino-Veltman 

In comparison to the above program the Mainz method may be characterized by the 
following list of statements 

• integrate in a special frame (usually the rest frame of one of the external particles) 



4 



• no Feynman parameters 

• do the integrations in Minkowski space 

• tensor integrals are worked out by direct integration 

The traditional integration methods become progressively more intractable when more 
and more mass scales are introduced into a problem. This has led many authors to intro- 
duce quite radical approximations such as setting the ratios m^fmt and m z /m t to zero in 
their loop calculations without being able to say much about the error invoked by such an 
approximation. Contrary to this the Mainz integration methods are ideally suited for the 
most general mass case with as many different massive loop particles as possible. In fact 
Standard Model calculations with m t ^ m-w 7^ vnz 7^ Trig are quite welcome. In the other 
extreme one generally has no problems taking the mass zero limit in a loop contribution if 
the diagram does not contain any infrared (IR) or mass (M) singularities. In the presence 
of IR/M singularities the corresponding integrals have to be regularized by subtraction as 
will be discussed later on. 

In order to provide a measure of the scope of what XL OOPS can do at the moment 
and what XLOOPS hopes to be able to do in the future we have drawn a representative set 
of loop diagrams in Fig|J which are labelled as: Al (two-point one- loop), A2 (two-point 
two-loop), Bl (three-point one-loop), B2 (three-point two-loop), CI (four-point one-loop) 
and C2 (four-point two- loop). In Table 1 we list the status of these loop calculations as 
concerns their availability in the present version XLOOPS 1.0 and, if not implemented in 
XLOOPS 1.0, we comment on their stage of development. 



loop class 


status and availability 


2-point 1-loop 


XLOOPS 1.0 


2-point 2-loop 


XLOOPS 1.0 


3-point 1-loop 


XLOOPS 1.0 


3-point 2-loop 


UV divergences : ready 

IR/M divergences : starting 

tensor structure : conceptually ready 


4-point 1-loop 


ready for implementation 


4-point 2-loop 


starting 



Table 1: Present status of one- and two- loop calculations in XLOOPS 

3 Parallel and orthogonal momentum space 

A very important technical concept in our loop calculations is the splitting of total d- 
dimensional space into a parallel space spanned by the external momenta and an orthogonal 
space spanned by its orthogonal complement. Symbolically we write 
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A1 ~~ C^^ - 


A2 








B1 




B2 






C1 


C2 





















Figure 3: Representative set of diagrams to be treated in XLOOPS 

total space = parallel(||) space + orthogonal(_L) space 

The orthogonal space components of the loop momenta have no orientation in the real 
world defined by parallel space. For scalar integrals this implies a spherical symmetry in 
orthogonal space which can be nicely exploited. Tensor integrals involving orthogonal com- 
ponents can easily be evaluated since their expansion along resultant tensors involves only 
combinations of orthogonal space metric tensors. These can be easily processed further. 
Also, by splitting 7-matrices into their orthogonal and parallel space components almost 
all the tedious 7-matrix algebra in fermionic amplitudes can be avoided. 

Let us illustrate the advantages of working in orthogonal space with the help of a few 
examples. We begin with the two-point function with any number of loops as drawn in 
symbolic fashion in Fig.|j. 




Figure 4: Two-point function with any number of loops 

The external momentum and thus parallel space is defined by the momentum q^. The 
dimensionality of the two respective spaces is thus given by 

parallel space (||) : 1-dimensional 
orthogonal space (_L) : (d — l)-dimensional 
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The metric tensor splits up into a parallel piece g\\ and a orthogonal piece g^ v 
according to = g\\ + g^ u . The respective components of the metric tensor can easily 
be constructed, viz. 

i = QnQv/q 2 (1) 
gi„ = 9^-q^/q 2 (2) 

Quite naturally the two components of the metric tensor can be used to project out the 
parallel and orthogonal components of any d- dimensional vector. 

In a similar way one can split up the 7-matrix into its parallel and orthogonal compo- 
nents according to 

7, = 7ji + 7^ (3) 

By using the metric projectors it is easy to see that one has t| = iq^/q 2 and 7^ = 
In ~ iq^/q 2 - Note also that 7JI and 7^ anticommute, i.e. one has 



7/i) 1v 



(4) 



The anticommutativity of the two respective 7-matrix components results in an extremely 
expedient way of dealing with the 7-matrix algebra when evaluating fermionic Feynman 
diagrams. We shall return to this point later on. 

In the case of the n-loop three-point function one has two independent external mo- 
menta q and p that define parallel space. A symbolic representation of the n-loop three- 
point function is drawn in Fig|| 




Figure 5: Three-point function with any number of loops 

Assuming that the momentum q is time-like one chooses the rest frame of q as the 
reference frame (if q is space-like one chooses the rest frame of one of the other momenta). 
In this frame the momenta q and p are represented by q = (q°, 0, 0, 0) and p = (po,p z , 0, 0) 
such that the three-momentum p defines the ^-direction in parallel space. The dimensions 
of the two respective spaces is now given by 

parallel space (||) : 2-dimensional 
orthogonal space (_L) : (d — 2)-dimensional 

As before the metric tensor splits into the the two components 

9i*v = A + g» v ( 5 ) 
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where one has 

II 1 1 ( ^ 2 \ ( R \ 

3% = -^W» - ^ \PuPu - —^Wv J ( 6 ) 

and gj[ v is given by the complement g^ v — gj w . 

We briefly side-track to remind the reader that the construction of a parallel and an 
orthogonal space is a familiar concept used also in other areas of particle physics even 
if different names are attached to the two respective spaces. For example, one talks of 
the longitudinal and transverse component of a four-vector current with respect to its 
momentum q. The transverse component is then defined through the relation q^J^ = 0. 
In a similar vein in HQET one chooses a time-like four-velocity with v 2 = 1 which serves 
to define longitudinal and transverse components of a general four-vector a M according to 
aj[ = a • and aj^ = a M — a • vv^. In the dynamical formulation of HQET this splitting 
up of four- dimensional space-time becomes an important dynamical concept in as much 
as the covariant derivative is split up into its longitudinal piece Z?j| = D ■ vv^ and its 
transverse piece = D^ — D-vv^. In the v rest frame = (1, 0, 0, 0) the longitudinal and 
transvere covariant derivatives reduce to the time-derivative Dq and the three-derivative 
D, respectively. The I/tuq expansion of the theory must then be organized in such a way 
that in the rest frame of v higher order derivatives involve only the three-derivative D 
whereas the time-derivative D must appear only to first order [ |TU| . 

In the two examples discussed before (two- and three-point function) orthogonal space 
was (d — 1)- and (d — 2)-dimensional. The continuation of space-time away from d = 4 
facilitates the regularization of ultraviolet (UV) and IR/M singularities. All one-loop 
calculations in XLOOPS are in fact carried out in d ^ 4 dimensions in order to be able to 
control the UV and IR/M divergences. As it turns out the two-loop integrations cannot 
be done in d ^ 4 dimensions in the generality that is needed for the implementation in 
XLOOPS. Instead the divergent pieces of the integrands are first identified. One then 
subtracts and adds auxiliary integrand functions that have the same degree of divergence 
but which are sufficiently simple to allow for their integration in d ^ 4 dimensions. The 
divergent behaviour of the original integrand is thereby eradicated such that its numerical 
integration can be performed in d — 4 dimensions. The divergent auxiliary function is 
then integrated using dimensional regularisation to the required accuracy in powers of 
£ = (d-4)/2. 

In the UV case the auxiliary integrand functions are constructed from the same Feyn- 
man diagram which the auxiliary integrand functions are sought to reproduce in their UV 
behaviour, however, with some of the momenta and masses in the diagram set to zero. It 
is clear that one has to be careful to avoid introducing spurious IR/M divergences in the 
course of this procedure. The construction of the requisite auxiliary integrand functions 
has been fully automated and is implemented in XLOOPS. 

If a diagram contains IR/M divergences one proceeds in a similar fashion. However, 
the construction of the requisite auxiliary integrand functions is not as simple as in the UV 
case. As of yet we have not been able to solve the problem of constructing the auxiliary 
integrand functions in all generality. There is a certain amount of progress, though, in 
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that we have been able to demonstrate the feasibility of this approach for some specific 
and important cases. The two cases treated so far are the two- loop 0{a 2 s ) contributions to 
the b — > c transition at zero recoil and the two-loop contributions to Z — > bb. 

There are altogether 13 two-loop contributions to the transition b —>■ c. According to 
the Lee-Nauenberg theorem the sum of the 13 two-loop contributions are IR finite at zero 
recoil because the accompanying Bremsstrahlung diagrams are kinematically suppressed 
at this point and can therefore not be called upon for IR cancellation. But taken sep- 
arately single diagrams generally do have IR singularities which have to be regularized. 
The requisite set of auxiliary integrand functions have been explicitly constructed. After 
subtraction the remaining IR finite contributions were integrated analytically. We mention 
that the corresponding one-loop calculation had been done in the early 80's for any value 
of the momentum transfer variable |TTJ. The first complete analytical zero recoil two-loop 



calculation was published this year by Czarnecki and Melnikov[12, 13 




Figure 6: Parallel ladder two-loop diagram contributing to Z — > bb. 

The second example concerns the determination of the 0(g 2 a s ) mixed electroweak and 
QCD two-loop contributions to Z — > bb. These radiative corrections are important for high 
precision tests of the Standard Model. In Fig.|6| we show one of the contributing diagrams 
which is of "parallel ladder" type. In our first exploratory calculation we neglected spin 
effects and calculated only the scalar parallel ladder diagram. The diagram is IR divergent, 
and, for rrif, = 0, also M divergent. The IR/M behaviour of this diagram and a few other 
parallel ladder diagrams was improved in the manner described above. The remaining 
IR/M finite integration was computed numerically using the two-fold integral representa- 
tion of massive two-loop contributions which is also available in XLOOPS. The auxiliary 
integrand functions were then integrated in d ^ 4 dimensions obtaining its IR/M divergent 
piece as well as its finite contribution. In order to check on the results of this numeri- 
cal evaluation the contribution of the same diagram was then calculated using completely 
different integration methods, namely the improved small momentum expansion method 



developed in Refs.JTj]. The numerical results were compared and satisfactory agreement 
was found [ |I~5|j . 

For the remaining part of this presentation we presume that the UV and IR/M im- 
provement program for the integrands has been carried through such that there is no need 
for regularization. The remaining discussion will be centered around integration methods 
in d = 4 dimensions. 
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4 Tensor integrals and strings of 7-matrices 



In order to illustrate the advantages of using parallel/orthogonal space techniques we return 
to the two-loop two-point case. In Fig.[7]we have drawn the generic diagram of the so called 




Figure 7: Two-loop two-point function: the master topology 

master topology which is the most difficult diagram in the two-loop two-point category. 
Included are also the flow of the loop momenta I and k. The propagator line 3 contains 
both loop momenta I and k and will therefore be referred to as mixed propagator line. 
The other propagator lines 1, 2 and 4, 5 contain only one loop momentum, namely / and 
k, respectively. Assuming that q is time-like, the diagram is evaluated in the g-rest frame 
where g M = (q°, 0,0,0). (If q is space-like the corresponding contribution can be obtained 
by analytic continuation from the time-like region.) The loop momenta split up into their 
parallel and orthogonal space components according to / M = (7°, l±) and k^ = (k°, k±). 

Assume now that there is a 7-matrix string associated with the diagram. Remembering 
that parallel and orthogonal 7-matrices mutually anticommute and that they anticommute 
with 75 the string can always be arranged into pieces of the form 

lTlo°t---t (7) 

where the 77,5 7s's have been moved to the left end of the string and the no 70's to second 
to left. The products of 7s's and 70's can then be simplified by using 7I = 1 and 7q = 1. 
If one has an even string of ^-matrices the string of 7 -matrices can be replaced by 
products of orthogonal space metric tensors whose weights are determined by the trace of 
the 7 -matrices. For those -y^-matrices that are contracted with the orthogonal space loop 
momenta /j_ and k± an appropriate symmetrization has to be carried out w.r.t. multiply 
occurring loop momentum indices. That is, returning to the covariant description, one has 
the replacement 

it ■ ■ ■ it • • • + ■ ■ • + symmetrization (8) 

Odd strings of ^-matrices can be seen to give zero results after integration. The case that 
the two-loop two-point master diagram carries external tensor indices will not be discussed 
further here but similar simplifications on the 7-matrix algebra can be seen to occur. 

Returning to Eq.(^) the Dirac-string object can be replaced by Zo^o — l± • k±. The 
orthogonal space scalar product l± ■ k± contains the cosine of the polar angle between l± 
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and k±. For further processing it is advantageous to get rid of the polar angle dependence 
through cancellation with the mixed propagator P 3 = (I + k) 2 — mf . This can be achieved 
by writing 

l± • k± = \ (ll + k 2 -P ± -kl-ml + 2l k - P 3 ) (9) 

The mixed propagator factor P3 cancels against the corresponding denominator factor and 
one remains with a factorizable two-loop contribution. The other terms are genuine two- 
loop contributions. However, they no longer possess a cos# dependence in the integrand 
numerator. They can therefore be further processed as in the scalar case to be discussed 
in the next section. 

If there are external indices attached to a diagram one has to solve the problem of how 
the loop contributions distribute themselves onto the various covariants that determine 
the covariance structure of the diagram. Again this problem is solved very elegantly using 
parallel/orthogonal space techniques. As an illustration we choose the simple example of 
the two-point function carrying two external Lorentz indices corresponding to the self- 
energy of an off-shell boson. The covariance structure of the resulting contribution has the 
form Cig^ u + C^q^qu- The invariant C\ can be picked out by calculating the orthogonal 
component of the two- loop integral since q^q v has no orthogonal space component. After 
having calculated C\ the invariant Ci can then be obtained by considering the parallel 
space component of the two-loop integral. In the language of contractions the invariant 
C\ is obtained by contraction with since q^q v g^ = 0. Once C\ is known the invariant 
C2 can be obtained by contraction with the parallel space tensor q^q u . Admittedly the 
example that we have discussed is of a rather simple nature but it still serves to illustrate 
the important point that the parallel/orthogonal space method allows one to efficiently 
evaluate tensor loop integrals. Much of the tedious algebra of matrix inversions needed 
in the Passarino-Veltman method of calculating tensor integrals is avoided since the use 
of parallel/orthogonal space techniques lead to an almost complete diagonalization of the 
problem from the outset. 

5 The two-loop two-point function: 
the master diagram 

Here we give a brief review of the derivation of Kreimer's two-dimensional integral rep- 
resentation m for the scalar master diagram (see Fig. [7]). This is mainly intended to be 
a warming-up exercise before we discuss the non-planar three-point function in the next 
section. 

The master diagram depends on one external momentum, q, and five masses m, n (i = 
1,...,5). In the scalar case, it is UV-finite, and can therefore be calculated in d — 4 
without any subtractions: 

/ = J(g 2 ;m 1 ,...,m 5 ) = / / d 4 k d 4 / - - i - - , (10) 
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where 



P 1 = I 2 - m \ + ip 

P 2 = (I - qf — ml + ip 

P 3 = (k + 1) 2 - ml + ip 

P 4 = k 2 — m\ + ip 

P$ = (k + q) 2 - ml + ip 



For the time being we shall assume that the external momentum is time-like (q 2 > 0), 
so that we can work in its rest frame, in which q = (qo,0), k = (fa,k±) and I = (Iq,1±). 



Introducing kj_ 
propagators as 



\k±\, l± = \l±\ and the polar angle 9 between kj_ and l±, we rewrite the 



2 7 2 

uJ 3 -k ± 



Pi 

p 2 



I 2 



— 1 1 — 2k±l± cos 9, 



P 4 
P 5 



2 i 2 

u 4 - k ± , 

2 i 2 
^5 - k ±i 



<JJ2 
UJ 5 



\-m\ + ip 



\j{k~ qo) 2 -ml + ip 



(fa + 1 ) 2 -ml + ip 



k 2 -m\ + ip 



(fa + Qo) 2 -ml + ip. 



:i2i 



The important points to note here are that only the mixed propagator, P 3 , depends on 9, 
and that the differences Pi — P2 and P4 — P5 only depend on the parallel space variables 
fa and Iq. After performing three trivial angular integrations, the integral fllPf) becomes 



x 



7T~ 
1 

-1 



dfa 
dcos# 



dl 







1 1 1 

+ 



k 2 6k \ 



11 dZ 



P1P3P4 P1P3P5 P2P3P5 P2P3P4 



(13) 



Let us concentrate on the first termQ in the parantheses on the second line of Eq. (|T3| ) . 
The integral over the orthogonal space variables can be regarded as a vacuum integral with 
"masses" u>i, 0J3, uj± in d = 3. The z integral is elementary: 



k\ dfcj 
dfa 
dfc 



I 2 , dl 



Z1 dZ, 



P1P4 
1 



PxP 4 2k 1 1 



k\l\ 



1 1 
dcos#— 
1 P3 

log 



(f j, - k ± ) 2 - u\ 
(l ± + k ± ) 2 -ul 



oo d/± (/i-^)(A;i-^ 4 2 ) 



log [l ± + k ± + u 3 ] 



(14) 



1 We are cheating slightly because the integral of 
cancels if the four terms in (|l^) are taken together. 



by itself diverges. However, the divergence 
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Now, both the k± and l± integrations can be performed by closing the integration contours 
in the complex plane (taking care not to cross the cut of the logarithm) and applying 
Cauchy's theorem, which gives: 



(2ni) 



log [iUi + UJ 3 + Co> 4 ] 



(15) 



Collecting the four terms in flT3] ) and expressing everything in terms of the dimensionless 
variables x = lo/qo, V — ^o/<?o ; we finally obtain: 



47f 4 /-oo /-oo 

— / dx dy — 

q z J-oo J-oo w{ 



1 



w 2 w\ 



with 



log 



(wi + w 3 + w 4 )(w 2 + w 3 + w 5 ) 
(w 2 + w 3 + w 4 )(wi + w 3 + w 5 ) 



(16) 



Wi 



\ 



x 2 — 



ip 



w 3 



w 5 



\ 



{x + yf 



+ ip 



w 2 



W4 



\ 



(x-r 



+ ip 



\ 



2 m 4 , • 

y — 5- + w 



\ 



(y + l) 2 -^ + *p. 



(17) 



The representation ( jl6|) is valid for arbitrary masses, and it gives the correct analytic 
continuation to space-like external momenta. Except for a few very special threshold 
configurations, it is possible to take the limit p — >• under the integral sign, yielding 
an integrand which is continuous in the entire (x, y) plane and can easily be integrated 
numerically. An attractive feature of this result is that we do not have to distinguish 
whether q 2 is below, in between, or above the thresholds: one simple formula covers all 
cases. Eq.(|TBp also gives the correct imaginary part. 

Tensor integrals (of rank greater than one) require subtractions to make them UV- 
convergent. XLOOPS automatically performs the necessary subtractions || and then 
generates two-dimensional numerical integrands for the finite remainders that are very 
similar to the one derived here for the scalar case. The only difference is that they contain 
additional factors of x, y, or Wi in the numerator which are not difficult to deal with. 



6 The two-loop three-point function: 
the crossed ladder configuration 

Starting from the two-point master topology, one obtains the crossed ladder topology by 
attaching the third leg to the central propagator, as shown in Fig.|[ This is the most 
complicated two-loop three-point topology because, no matter how the loop momenta k 
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Figure 8: The two-loop crossed vertex diagram 



and I are routed, one cannot avoid having two mixed propagators that depend on the scalar 
product k ■ I. Another way of drawing this topology is shown in Fig.|TT|. 

In this section, we will give a brief and somewhat simplified account of how a two-fold 
integral representation for this diagram can be obtained. For a more detailed explanation, 
and for other two-loop three-point topologies, we refer to [[16|, ||. 

The scalar integral is defined by 

V = V(qlqlq 2 ;m 1 ,...,m 6 ) = f f d 4 k d 4 Z 1 , (18) 

where q = q\ + q 2 and the propagators are 

Pi = (k + I - qxf - m\ + ip 



P 2 = (k + l + g 2 ) 2 — m\ + ip 

P3 = {I - qif -ml + ip 

Pa = (k + q2) 2 - ml + ip 

P5 = I 2 — ml + ip 

P 6 = k 2 -ml + ip. (19) 

We use a frame where q = (ei+e 2 ,0,0), q x = (ei,g 2 ,0), q 2 = (e 2 ,-g z ,0), k = {k 0l ki,ki_) 
and / = (lo,li,l±). This time, k± and l± are two-component vectors. Introducing polar 
coordinates in the orthogonal space and carrying out one trivial angular integration, we 
find that 

/d'*/<rt... 

/oo roo poo /"CO roo roo /■27r 

dk / d/ / dh / dh / k± dk ± / l ± dl ± d6 . . . 

-00 J-00 J~oc J-00 JO JO Jo 

00 _ /"oo _ roc roo roo roo rliv 

dk / dZ / ds dt dh dh d9 . . . , (20) 
2, J— 00 J —00 Jo Jo J —00 J— 00 Jo 



7T 



where s = l\ and t — k\. In the last line, we have replaced ko and Iq with the shifted 
variables ko = k — ki and Iq = Iq — l±. As a result, all the propagators become linear in ki 
and li, e.g. 

Pq = kl — k\ — t — ml + ip = k 2 , + 2fc ^i — t — ml + ip . (21) 
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Thus, the k\ and l\ integrations can be performed by closing the contours in either the 
upper or lower half plane and applying Cauchy's theorem. In order to understand the 
general structure of the result of using Cauchy's theorem, let us decompose the integral 
into (2x2x2 = 8) partial fractions of the form 

dJfei f°° d/ i n v7 Kn i V = ^ ■ ( 22 ) 

./-oo (h - - z 2 )(k 1 + h + z 3 ) z 1 + z 2 + z 3 

This equation holds when the signs of the imaginary parts of z\, z 2 and z% are all equal. 
Otherwise, the integral on the left hand side vanishes. Now the imaginary parts of the 
Zj's all come from the infinitesimal ip, divided by linear combinations of ko, Iq and the 
external momenta (cf. eg. eq. fl2lD ). As a consequence, each term of the partial fraction 
decomposition only gives a non-vanishing contribution inside a finite, triangular region in 
the (ko, Iq) plane. The union of those triangles is shown in Fig. |9|. 

After a further partial fraction decomposition wrt. cos 9, the ^-integration can be carried 
out using 

r2n I 27T 

dO- A ^ ~ n = , (23) 

A-Bcosd ^A 2 - B 2 

(and being careful to pick the correct branch of the square root on the right hand side!). 

Having performed the k±, l\ and 9 integrations^, one obtains the following intermediate 
expression: 

n ft ^ i 

V = E / = IT • (24) 

(djt + bj + Cjs) 2 - Ast i= i (aijt + b i;j + djs) 



The quantities Cj, a>j, bj, Cj, 5^, bij and c^- are all rational functions of ko, lo, q% and rn.j. 9j 
is a product of step functions restricting k and Zo to one of the triangles discussed above. 

By means of Euler's change of variables, it is possible to get rid of the square root in ( pip 
and perform both the s and t integrals analytically, leading, eventually, to dilogarithms and 
Clausen functions. Although the final expression is too long to be listed here in explicit 
form, it is a sufficiently well-behaved function to allow the remaining ko and Iq integrations 
to be evaluated numerically without problems. 



7 Checks on two-loop three-point functions 

We have performed a number of checks on the integration methods that go into XLOOPS 
as well as on the performance of XLOOPS itself. Some of these checks have been discussed 
before, as e.g. checks on the two- loop 0(a^) contributions to the b —>■ c transition at zero 
recoil. 

As concerns the two-loop three-point function the checks divide up into checks on the 
correctness of the theoretical structure of the results and numerical cross-checks against 
the results of other authors for specific mass cases. For example, we have checked that the 
two-fold integral representations of the parallel and crossed ladder contributions have the 

2 For technical reasons, the 9 integration in is done first. The result is, of course, the same. 
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Figure 9: Integration region in the (k 0} l ) plane 



• correct cut structure in the analytic plane as prescribed by the Landau- Cutkosky 
cutting rules 

• possess the correct anomalous threshold behaviour at the prescribed anomalous 
thresholds 

We have looked at all two-particle cuts and, with particular care, at all three-particle 
cuts and have found that the analyticity structure of our parallel and crossed ladder two- 
loop results are in complete agreement with the Landau-Cutkosky cutting rules. 

As a check on the correct anomalous threshold behaviour we ran our numerical program 
across an anomalous threshold. We start with the one-loop vertex contribution which 
possesses anomalous thresholds at 

1 + 2^1/^3 - lA - ~ f4, — (25) 

where (i = 1, 2, 3) 

= (raj + ml + mj - mf - p?)mj ^ 
2mim2m3 

As sample values we took the following set of masses and momenta (in arbitrary units) 

mi = m,2 = rris = 1 (27) 

p\ = 6, p 2 2 = 5, p 2 3 = -4 - Vl5 « -7.873 (28) 

In Fig.pl] we show the values of the imaginary and the real partj^of the vertex contri- 
bution in the vicinity of the anomalous threshold at p\ ~ —7.873. The plot clearly shows 
the logarithmic type singularity of an anomalous threshold at X — Xq which reads 

i ln(x — x ) + C(x) = i In \x — x \ ± tc9(x — x ) + C(x) (29) 



3 Because our vertex function contains an extra factor of i the role of the imaginary and real part is 
interchanged relative to what one is used to 
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Figure 10: Plot of real and imaginary parts of a vertex function in the vicinity of an 
anomalous threshold, one-loop (left) and two-loop (right) 



where C(x) is a complex- valued function that is regular in the vicinity of x = xq. We also 
checked on the anomalous threshold behaviour of the corresponding parallel ladder two-loop 
vertex correction. The two-loop vertex correction has an anomalous threshold at the same 
position p\ ?a —7.873 for any values of the masses 7714,1715, rag. The numerical evaluation 
of the two-loop graph showed the same behaviour around the anomalous threshold as the 
one-loop graph (see Fig.|H]). The numerical stability in the direct vicinity of the anomalous 
threshold was found to be satisfactory in both cases but was somewhat decreased in the 
two-loop case relative to the one-loop case. 

Further we ran our program for specific mass configurations and compared the numerical 
results with the results of other groups. For example, we evaluated the crossed ladder 
vertex graph for external momenta with q 2 = p\ = p\ 7^ and let the internal masses 
mi = mii = 1, . . . , 6) go to zero. When m = Ussyukina and Davydychev have shown 
that the crossed ladder configuration can be calculated analytically, and that the result 
turns out to be equal to the square of a one-loop graph jl7 . Our numerical m — ► results 
agree with the Ussyukina-Davydychev limit. It came as a pleasant surprise that numerical 
stability was retained in this limit. 

Next we took the specific momentum configuration p\ = p\ = and set all six internal 
masses in the crossed ladder configuration to m = 150 GeV. We compared our results with 
the corresponding calculation by Tarasov [18], who used an improved small momentum 



expansion using conformal mapping and Pade approximants. There was full agreement on 
the real and imaginary parts of the crossed ladder configuration above, below and close to 
threshold q 2 = Am 2 . In particular there was no noticeable decline in numerical stability 
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close to threshold which had been a troublesome region in the low momentum expansion 
method before it was improved by conformal mapping. It is not so simple to compare 
the time performance of the two approaches since the low momentum expansion method 
requires a great amount of preparatory work in that the coefficients of the low momentum 
expansion have to be calculated first. 

We also compared our results to an evaluation of the crossed ladder contribution to 
Z (off-shell) — > ti done by the Japanese KEK group Fujimoto et al.[19|. The relevant 



mass configuration is m x = m 2 = m 3 = m 4 = J q\ = Jq\ = 150 GeV (the would- 




Figure 11: Crossed ladder two-loop contribution to Z (off-shell) — > tt. 

be mass of the top quark in 1994) and = rriQ = 91 GeV. The contribution clearly 
corresponds to the off-shell Zti vertex with internal Z and t exchange. Again we found 
full numerical agreement for both the real and imaginary parts. The Japanese group used 
Feynman parametrization and reduced the number of necessary numerical loop integrations 
to three by analytic integration. For comparable accuracy the necessary computing time 
is considerably smaller in the Mainz approach. 

As a final example of the reliability of the massive two-loop crossed and parallel ladder 
integration techniques and the numerical integration techniques based on them we discuss 
a recent application of these methods to the calculation of the dominant two-loop correc- 
tions, of order 0(G 2 F Mfi), to the partial widths of the decays of a heavy Higgs to pairs 
of W- and Z-bosons, carried out in collaboration with B.A. Kniehl and K. Riesselmann 
20| . The calculation was done again using the Equivalence Theorem, i.e. the dominant 



longitudinal degrees of freedom of the gauge bosons were represented by the massless and 
spinless Goldstone boson excitations of the effective theory. Of the many contributing dia- 
grams we show only two representative diagrams of the class of crossed and parallel ladder 
configurations which are the most difficult to calculate. 

In the numerical evaluation of these diagrams we started off with small finite Goldstone 
masses which were then set to smaller and smaller values and finally to zero. We found that 
the numerical results were stable in this limit. After compounding the contributions of the 
many contributing diagrams and carrying out the renormalization in the on-shell scheme, 
the correction factor was found to be 1 + 14.6%(7nff/TeV) 2 + 16.9%(m^/TeV) 4 , indicating 
that for mn > ITeV the Standard Model ceases to be weakly interacting. Numerically, this 
result agrees with that of Ghinculov, who had used entirely different integration techniques 
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Figure 12: Two- loop parallel and crossed ladder contributions to H — ► ZZ in the effective 
theory. Solid lines: massive Higgs particles; dashed lines: massless gauge bosons. 



to evaluate the same decay process pl| . From the fact that the calculation involves a 
large number of massive two-loop diagrams that contribute to the process we believe that 
the numerical agreement of the two calculations provides for another strong test of the 
reliability of our integration techniques and the numerical routines based on them as well 
as those used by Ghinculov. We mention that a two-loop calculation of the fermionic 
contributions to H — > ZZ, WW is under way. 



8 Some sample calculations with XLOOPS 

This summer saw the premiere of a live performance of XLOOPS during a talk given by 
J. Franzkowski at the Klausurtagung of the Mainz Graduiertenkolleg "Teilchenphysik bei 
hohen und mittleren Energien" . For the performance he used a laptop and a projectable 
screen display. Nothing can beat the excitement of a live performance of XLOOPS either 
session or during a talk using a projector screen. We shall nevertheless try to 
capture some of the spirit of a live demonstration of XLOOPS by presenting two sample 
calculations done by XLOOPS. In the first example we calculate the on-shell one- loop 
self-energy of the neutral gauge boson Z resulting from the contribution of the top quark 
loop. In the second example we calculate one of the mixed 0(g 2 a s ) QCD/electroweak 
contributions to the two-loop flavour- changing self-energy of the s — > d quarks. 

At the beginning of an XLOOPS session one has to specify the model whose Feynman 
rules are to be used in the loop evaluation. In both of the examples to be discussed here 
one clicks on "Standard Model" in the "Model" menu which would call up the Feynman 
rules of the Standard Model. Also the user has to specify the numbers m and n for his 
choice of the m-point n-loop function. In the case of the one-loop self-energy of the Z 
he would then click on "two" external lines (m = 2) and on "one" loop (n = 1). The 
topology bar would then show the two possible two-point one-loop topologies. These can 



be discerned at the top of Fig.[L3|. One then clicks on the desired topology which would be 



the second of the two topologies in FigO. After having selected the topology the desired 



one- loop graph appears on the screen in an enlarged version as shown in Fig.|13. 

One then has to attach arrows to the lines in the graph according to the quantum 
number flow (wrongly pointed arrows would lead to an error message) and specify the 



19 



File Model Loops Ext. Lines Options 



Help 



Q O 



Evaluate Graph 



Subtopologies 



Export 




Evaluate E. Full 



E. More aose 



1 . ««««««« « « « « « « « « « « « >x 
K3WElci f Dirk Kreimer 
verEitiit Mainz , Gemnany 



MapleVR3/laopE .ma' ) : 



Figure 13: XLOOPS display of a one-loop Feynman diagram. 



particles in the graph by inserting the particle names in the particle entry fields. In the 
next step one has to specify the mass and momentum parameters. There is an "Insert 
Particle Properties" option as a menu entry which automatically substitutes the mass 
values provided by the latest version of the Particle Data Group listings. Mass values can 
of course also be introduced by hand. In our example we have m z = 91.187 GeV and 
m t = 176 GeV. The one-loop graph is evaluated in the rest system of the on-shell Z. Thus 
one chooses qio = 91.187 GeV. The calculation of the graph is then executed by clicking on 
the "E. (Evaluate) Full" bar. In this one-loop example the numerical answer would appear 
within seconds in the following format: 

[ 

CI = [ 

-16 

146.22583756233107509 I, .26077710144029877423*10 - 1219.7690445131992210 I 
], 

C2 = [.00096061186286253340551 I, 

-20 

.18797591473406467665*10 - .0080036871858407297497 1], 
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CI G(nul, nu2) + C2 ql(nul) ql(nu2)] 



The numerical output is given in terms of the two invariants C\ and C2 in the covariant 
expansion Cig fiu + C2q^ t q u - Each of the invariants is displayed as a list [a, b] which represents 
a series a/e + b + 0(e). In the present case the coefficients a are purely imaginary while the 
6's have a tiny real part, that would vanish if the causal ip's in the propagators were put 
to zero exactly. This corresponds to the fact that in the present momentum configuration, 
below the tt-threshold, the self-energy £ (which differs from the diagram itself by a factor 
of i) is purely real. 




Figure 14: XLOOPS input panel for the two-loop master diagram. 

The second example is drawn from a recent determination of the 0(g 2 a s ) two-loop 
contribution to the flavour changing self-energy for s —>■ d. Of the contributing diagrams 
we take one of the diagrams with a master topology as drawn in Fig.|T3. One proceeds as 



was described in the first example only that one specifies to "two" in the menu "Loops" . 
One recognizes some of the relevant two-loop topologies at the top of the display depicted 
in Fig.|T4l. In the present example one clicks on the master topology which is the fourth 
two-loop diagram from the left in the topology bar which is covered up in the display 



Fig.O. The relevant Feynman diagram pops up in which the quantum number flows and 



particle names have to be inserted as shown in Fig.[Tj. Finally, after specifying the particle 



21 



masses and the momentum of the external particle one initiates the numerical evaluation 
of the diagram: 

-6 

[CI = [ - .68103432226158921830*10 I delta3(tmul , tmu2) , 

-5 

.42790442626661980279*10 delta3(tmul , tmu2) 

-5 

- .40233768346700174054*10 I delta3(tmul , tmu2) , 
[7.1, [ - 4.98962 - 2.43436 I, .0015744 + .00190137 I]], 
.000026216132920579719619 delta3(tmul, tmu2) 

-6 

+ .60792896861672700256*10 I delta3(tmul , tmu2)] , 

-6 

C2 = [ - .68103432226158921830*10 I delta3(tmul , tmu2) , 

-5 

.42790442626661980279*10 delta3(tmul , tmu2) 

-5 

- .40233768346700174054*10 I delta3(tmul , tmu2) , 
[°/„l, [ - 4.99119 - 2.43249 I, .00176832 + .00185428 I]], 
.000026216121704383664719 delta3(tmul , tmu2) 

-6 

+ .59969113686867746914*10 I delta3(tmul , tmu2)] , 
C3 = [.00015525369186016884013 I delta3(tmul , tmu2) , 

- .00051268542153584540755 delta3(tmul , tmu2) 

+ .00072319063843461738822 I delta3(tmul , tmu2) , 
[°/„l, [ - 845.404 - 954.737 I, .742148 + .926344 I]], 
3.5949563175956949478 delta3(tmul, tmu2) 

+ 96495.670995961263803 I delta3(tmul , tmu2)] , 
C4 = [.00015525369186016884013 I delta3(tmul , tmu2) , 

- .00051266860707202340491 delta3(tmul , tmu2) 
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+ .00072847640460045290512 I delta3 (tmul , tmu2) , 
[7.1, [ - 844.468 - 954.961 I, .576421 + .882752 I]], 
10.877778389998066059 delta3(tmul, tmu2) 

+ 103995.51572482192447 I delta3(tmul , tmu2)] , 
CI (1 &* ONE) + C2 (1 &* Dg5) + C3 (1 &* DgO) + C4 &*(1, Dg5, DgO)] 

-6 

7.1 := .85129290282698652282*10 delta3(tmul , tmu2) 

The results are given in terms of the four invariants C\, C*2, C3 and C4 in the rest 
frame decomposition C\uu + C^u^^u + C^u^qU + C4UJ5J0U. For each of the invariants 
Cj, the data are shown as a list [a, b, [c, [d, e]), /]. Their meaning is as follows: Cj = 
a/e 2 + b/e + cx (d±e) + f + 0(e). f is the piece of the finite part of Cj that is calculated 
analytically, d is the part that is obtained by numerical integration, and e is an estimate 
of the accuracy of the numerical integration. The overall colour factor of the diagram 
delta3(tmul , tmu2) is displayed separately. It is calculated from the rules of the SU (N) 
colour algebra which are known to XL OOPS. 
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Figure 15: List of four-point two-loop topologies. 
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As indicated in Fig.^| and Table 1 we have started work on massive four-point two- 
loop functions. As a preparatory step we have classified all the two-loop topologies and 



generated their diagrams. All the four-point two-loop topologies are shown in Fig. 15 



We shall refer to the p-th diagram in the q-th row by (p,q). The calculation of diagram 
(2, 4) was done first because of its implication for a two-loop chiral perturbation theory 
calculation of 7 + 7 — ► n + 7T. The calculation of the scalar diagram (2,4) was recently 
completed. The techniques of parallel and orthogonal space integrations used in the four- 
point case are similar to the three-point two-loop case only that parallel space is now 
3-dimensional. Six of the eight needed integrations were done analytically. The remaining 
two-fold integral representation was found to be suitable for numerical processing. The 
calculation of diagrams (3, 3) and (7, 3) is well under way. 



9 Present status of XLOOPS and outlook 

The presently available version of XLOOPS is a demonstration version which runs under 
the name of XLOOPS 0.9. It is meant to whet the appetite of the user. The demo version 
features all Standard Model Feynman rules including QCD. It allows one to calculate 
one-loop graphs up to and including three-point functions. It is available at the web-site 
http: / /wwwthep. physik.uni-mainz.de/ ^xloops| . 

A full version XLOOPS 1.0 has been completed and is currently being tested. It is 
accompanied by a 110 page manual which can be obtained at the above web-site. XLOOPS 

1.0 and the manual should be available in a few weeks time. There are some additional 
features as Feynman rules for the two Higgs doublet model. One of the future items of 
XLOOPS is the automatic generation of Feynman diagrams from a given set of Feynman 
rules. One would then like to avail of subsets of the Standard Model Feynman rules. This 
feature is already incorporated in XLOOPS 1.0 in that one can select the Feynman rules of 
QED and QCD. XLOOPS 1.0 incorporates first two-loop features in the form of two-loop 
two-point functions including their tensor decomposition. 

As time goes on more features will be added to XLOOPS. We plan to issue XLOOPS 

1.1 by the end of this year. It will contain a one-loop four-point function routine including 
tensor structure capabilitiies. Also XLOOPS 1.1 will be able to automatically generate 
all one- and two-loop Feynman diagrams that contribute to a given process including a 
drawing routine which generates Postscript output. Also we plan to install a plotting 
routine that will generate plots when input parameter values of masses or momenta are 
varied continuously over a range of the input parameter values. There are also plans to 
write a longer review on the theoretical background with a detailed description of the 
calculational techniques that go into XLOOPS . 

Work is going on on the two-loop three-point function which is ready conceptually. 
More Feynman rules will be added such as those of the minimal supersymmetric Standard 
Model and chiral perturbation theory. Further in the future lies the implementation of 
the two-loop four-point function which is presently under study. Also we plan to install 
automatic renormalization capabilities. The ground is prepared for renormalization in 
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that the one- loop routines calculate to 0(e) accuracy such that one-loop counter terms 
can be computed automatically. In the very far future we hope to provide the possibility 
to evaluate complete processes including phase space integration routines for tree graph 
contributions. All of this needs a great deal of man power. We hope that we can count on 
the support of the community in this enterprise. 
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